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The emission of nuclear clusters is investigated within the framework of isospin dependent lattice 
gas model and classical molecular dynamics model. It is found that the emission of individual 
cluster which is heavier than proton is almost Poissonian except near the transition temperature at 
which the system is leaving the liquid-vapor phase coexistence and the thermal scaling is observed 
by the linear Arrhenius plots which is made from the average multiplicity of each cluster versus the 
inverse of temperature in the liquid vapor phase coexistence. The slopes of the Arrhenius plots, 
i.e. the "emission barriers", are extracted as a function of the mass or charge number and fitted by 
the formula embodied with the contributions of the surface energy and Coulomb interaction. The 
good agreements are obtained in comparison with the data for low energy conditional barriers. In 
addition, the possible influences of the source size, Coulomb interaction and "freeze-out" density 
and related physical implications are discussed. 
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I. INTRODUCTION 



Fragmentation phenomenon is an important subject in 
many different research fields, such as in nuclear physics 
and cluster physics etc. It is found that fragmentation 
phenomenon has a close relation to liquid-to-gas phase 
transition in heavy ion collision physics 0, 0, 0, 0, Q 
and cluster physics [g, 0- In heavy ion collisions, the 
fragmentation can occur from a short-lived hot excited 
nuclei with moderate temperature which can be formed 
by the collisions between heavy ions at low-intermediate 
energy. Usually the hot nuclei finally de-excite by the dif- 
ferent decay modes, such as the light particle evaporation 
and the emission of multiple intermediate mass fragment 
(IMF), i.e. multifragmentation. Even though extensive 
studies on fragmentation have been carried out experi- 
mentally and theoretically, it is still an open question to 
clarify definitely whether the fragmentation is statistical 
or dynamical, sequential or simultaneous. Among such 
efforts, the reducibility and thermal scaling in multiple 
fragment emission process have been explored |g, 13, Il0| 
and seems to show one a possible interpretation picture 
to look and understand the multifragmentation. Origi- 
nally, Moretto et al. observed that the experimental Z- 
integrated fragment multiplicity distributions P™ are bi- 
nomially distributed, i.e. 



in each transverse energy (E t ) window, where n is the 
number of emitted fragments and m is interpreted as the 
number of times the system tries to emit a fragment. 
The probability of emitting n fragments can be reduced 
to a single-particle emission probability p which gives 
linear Arrhenius plots (i.e. excitation functions) when 
ln(l/p) is plotted vs \j \[W t . By assuming a linear re- 
lationship between \fE t and temperature T, Moretto et 
al. proposed that the linearity of the observed ln(l/p) 
vs \j\fEt can be interpreted as a thermal scaling of the 
multifragment process |a H, [Hj ■ In this case, these lin- 
ear Arrhenius plots suggest that p has the Boltzman form 
p ex e ' with a common fragment barrier B. However, 
since the binomial decomposition has been performed on 
the .Z-integrated multiplicities, typically associated with 
3 < Z < 20, the Arrhenius plot generated with the result- 
ing one fragment probability p is an average over a range 
of Z values. However, some debates and comments on 
the above binomial distribution and the thermal scaling 
were also raised [H 111 13 

Later, instead of analyzing for Z-integrated multiplici- 
ties, the behavior of individual fragment species of a given 
Z for higher resolution experimental data was investi- 
gated and found that the n-fragment multiplicities P(n) 
obey a nearly Poissonian distribution [151 ll(i| , 
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where n is the number of fragments of a given Z and 
the average value (n) is a function of the total transverse 
energy E t , and were thus reducible to a single- fragment 
probability proportional to the average value (n) for each 
Z. Similarly the (n) is found to be proportional to e~ B l T 
providing that T ex \fE~ t , i.e. there exists also a ther- 



mal scaling law. More recently, the common features 
of Poissonian reducibility and thermal scaling can also 
be revealed in percolation and the Fisher droplet model 

0m. 

In the present work, we would like to make a theo- 
retical re-examination on the Poissonian reducibility and 
its thermal scaling rather than the bionomial reducibility 
and its thermal scaling. Unlike experiment, we will adopt 
the true temperature as model input to check the Poisso- 
nian reducibility and thermal scaling in the frameworks of 
the isospin dependent lattice gas model (I-LGM) 19, 20] 
and followed by the isospin dependent classical molecu- 
lar dynamics (I-CMD) of Das Gupta and Pan [2l|. By 
investigating the variances and average multiplicities of 
cluster multiplicity distributions as a function of temper- 
ature, we will illustrate that the Poissonian reducibility 
and its thermal scaling is valid for the fragment emis- 
sion in the low temperature side in the framework of the 
above thermal equilibrium models. 

The paper is organized as follows. Firstly, we introduce 
the models of I-LGM and I-CMD in Sec. II; In Sec.III, 
the results and discussions are presented. First we show 
some results to support that there exists a Poissonian 
reducibility in the cluster production away from the liq- 
uid gas phase transition by investigating the ratio of the 
dispersion of multiplicity distribution to its mean value 
and the Poissonian fit to the multiplicity distribution of 
individual clusters. Second we plot the Arrhenius-type 
plots and find the thermal scaling is valid only in the 
liquid vapor coexistence phase. Further we extract the 
"emission barrier" sorted by the different mass number, 
light isotope, and charge number in I-LGM and I-CMD 
and use the formula embodied with the contributions of 
the surface energy and Coulomb interaction to make a 
systematic fit. The dependence of the model, the source 
size, Coulomb interaction and "freeze-out" density are 
presented. Finally, a summary and outlook is given in 
Sec. IV. 



II. DESCRIPTION OF MODELS 

Originally, the lattice gas model was developed to de- 
scribe the liquid-gas phase transition for atomic system 
by Lee and Yang [22] • The same model has already been 
applied to nuclear physics for isospin symmetrical sys- 
tems in the grand canonical ensemble l23| with a sam- 
phngof the canonical ensemble [H HaM IH |H 
1271 128 . l29l l30| . and also for isospin asymmetrical nu- 
clear matter in the mean field approximation [3l|. In 
this work, we will adopt the latti ce g as model which was 
developed by Das Gupta et al. j^J,|2£j. I n addition, a 
classical molecular dynamical model of Das Gupta et al. 
|2ll ] is also used to compare with the results of lattice gas 
model. For completeness of the paper, here we make a 
brief description for the models. 

In the lattice gas model, A {= N + Z) nucleons with 
an occupation number Si which is defined Si — 1 (-1) for 



a proton (neutron) or si = for a vacancy, are placed 
on the L sites of lattice. Nucleons in the nearest neigh- 
boring sites have interaction with an energy e SiS .. The 
Hamiltonian is written as 
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where Pi is the momentum of the nucleon and m is its 
mass. The interaction constant e SiS . is chosen to be 
isospin dependent and be fixed to reproduce the bind- 
ing energy of the nuclei J2l| : 
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Three-dimension cubic lattice with L sites is used which 
results in pf — -tPq of an assumed "freeze-out" den- 
sity of disassembling system, in which po is the normal 
nuclear density. The disassembly of the system is to 
be calculated at pf, beyond which nucleons are too far 
apart to interact. Nucleons are put into lattice by Monte 
Carlo Metropolis sampling. Once the nucleons have been 
placed we also ascribe to each of them a momentum by 
Monte Carlo samplings of Maxwell-Boltzmann distribu- 
tion. 

Once this is done the I-LGM immediately gives the 
cluster distribution using the rule that two nucleons are 
part of the same cluster if 



P^/2p- es iSi SiSj < 0, 



(5) 



where P r is the relative momentum of two nucleons and 
pi is their reduced mass. This prescription is evidenced 
to be similar to the Coniglio-Klein's prescription l32j| in 
condensed matter physics and be valid in I-LGM |19L |20L 
1261 1271 ] . To calculate clusters using I-CMD we propagate 
the particles from the initial configuration for a long time 
under the influence of the chosen force. The form of the 
force is chosen to be also isospin dependent in order to 
compare with the results of I-LGM. The potential for 
unlike nucleons is 20, 2lLl33l 
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In the above, r = 1.842/m is the distance between the 
centers of two adjacent cubes. The parameters of the po- 
tentials are p = 2, q = 1, a = 1.3, B = 0.924, and A — 
1966 MeV. With these parameters the potential is min- 
imum at ro with the value -5.33 MeV, is zero when the 
nucleons are more than 1.3ro apart and becomes stronger 
repulsive when r is significantly less than ro . The poten- 
tial for like nucleons is written as 



v PP (r)(r < r ) 
v P p(r)(r > r ) 



v np (r) -v np (r ), 

0. (7) 



This means there is a repulsive core which goes to zero 
at r and is zero afterwards. It is consistent with the fact 
that we do not put two like nucleons in the same cube. 
Essentially the classical molecular dynamics is based on 
the idea that the equation of state of classical particles 
interacting through attractive and repulsive Yukawa or 
Lennard-Jones potentials leads to an equation of state 
similar to a Van de Waals one. This kind of approach 
was first proposed in Ref. 34] to treat the nuclear col- 
lision with the similar potential as used in this paper. 
Extensive studies using the similar CMD to explore the 
liquid gas phase transition and fluctuations were also per- 
formed in some works, such as Ref. |35ll3a . I37ll3q . 

The system evolves for a long time from the initial 
configuration obtained by the lattice gas model under the 
influence of the above potential. At asymptotic times the 
clusters are easily recognized. The cluster distribution 
and the quantities based on it in the two models can now 
be compared. In the case of proton-proton interactions, 
the Coulomb interaction can also be added separately 
and compared with the cases where the Coulomb effects 
are ignored. 



III. RESULTS AND DISCUSSIONS 

In this paper we choose the medium size nuclei 129 Xe 
as a main example to analyze the behavior of individ- 
ual fragment emission during nuclear disassembly with 
the helps of I-LGM and I-CMD. In addition, the systems 
with A sys = 80 (Z sys = 33) and 274 (Z sys = 114) are 
also studied to investigate the possible source size de- 
pendence. One part of this work has been reported in 
Ref. pi] for 129 Xe in the I-LGM calculation with a fixed 
"freeze-out" density. In this work, a lot of new calcula- 
tions are largely included. 

In both model calculations, the "freeze-out" density pf 
is mostly chosen to be about 0.38 po, since the experi- 
mental data can be best fitted by p/ between 0.3po and 
0.4,oo in the previous LGM calculations [jjj,|4£j, which 
corresponds to 7 3 cubic lattice is used for Xe, 6 3 for A sys 
= 80 and 9 3 for A sys = 274. Under the condition of the 
fixed "freeze-out" density, the only input parameter of 
the models is the temperature T. In the I-LGM case, 
Pf can be thought as the " freeze-out" density but in the 
I-CMD case pf is, strictly speaking, not a "freeze-out" 
density but merely defines the starting point for time 
evolution. However since classical evolution of a many 
particle system is entirely deterministic, the initializa- 
tion does have in it all the information of the asymptotic 
cluster distribution, we will continue to call pf as the 
"freeze-out" density. 1000 events are simulated for each 
T which ensures enough statistics. 



A. Poissonian Reducibility 

One of the basic characters of the Poisson distribution 
Eq.(2) is the ratio cr 2 ./(rij) — > 1 where a 2 , is the vari- 
ance of the multiplicity distribution and (rii) is the mean 
value of the multiplicity distribution. The first step we 
are showing is this ratio. We obtain these ratios for clus- 
ters classified with different masses (A), light isotopes 
{ISO) and atomic numbers (Z) for the disassembly of 
129 Xe as a function of temperature in the framework of 
I-LGM and I-CMD with Coulomb in Fig. [l] Obviously, 
most of the ratios are close to one, which indicates that 
basically these cluster production obeys the Poisson dis- 
tributions, i.e. a cluster is formed independently from 
one another. Of course, we also notice that the values of 
protons are almost lower than the unique, i.e. it is nar- 
rower than the Poisson distribution. This could be due 
to protons can be easily separated from some unstable 
multi-nucleon clusters. In this case, this kind of proton 
production is obviously related to the parent cluster and 
then a narrower distribution of mixed protons could re- 
veal. On the contrary, some points are slightly larger 
than 1 and this behavior becomes obvious in the mediate 
temperature range, which is happening because at this 
temperature the system is leaving the liquid vapor phase 
coexistence [ljj, where the large fluctuation 41] makes 
the Poissonian reducibility broke-down. 

Besides the above Poisson condition is basically sus- 
tained, the excitation function of the average multiplic- 
ity of n-multiple individual cluster emission can be well 
fitted with the Poisson distribution. For some examples, 
Fig-Elshows the quality of the Poissonian fits to the aver- 
age multiplicity of n-multiple individual cluster emission 
in the different temperature for 129 Xe in the I-LGM case. 
In each panel of this figure, we first plot the probability 
of n- multiple emission species ( P(n) ) as a function of 
temperature (as shown by the open symbols), and then 
we connect the Poisson probabilities in the different tem- 
perature as lines in terms of Eq.(2) due to we know n 
and its average value (n) over all n-multiple emission in 
each temperature (as shown by the lines). Obviously, 
these Poissonian fits are quite good for almost species 
with Z > 2 over the entire range of T. The similar good 
Poissonian fit is overall obtained in the cases of I-CMD. 
As a consequence, the Poissonian reducibility is valid in 
the thermal-equilibrium lattice gas model or molecular 
dynamics, which illustrates that the cluster production 
is almost independent each other in the studied models. 



B. Thermal Scaling 

Naturally, we want to know if there exists a thermal 
scaling law in the thermal-equilibrium LGM and CMD 
models. To this end, the temperature dependence of 
the mean yield ((n)) of individual clusters is investi- 
gated. In order to compare a recent well-known thermal 
Arrhenius-type plot in nuclear multifragmentation phe- 
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FIG. 1: The ratio of (Ti 2 /(ni) for the clusters classified with mass, 
light isotope mass and atomic number as a function of temperature. 
The left panels are for the I-CMD calculation and the right for I- 
LGM with Coulomb. The symbols are illustrated on the figure. 



FIG. 2: The excitation functions of n-multiple cluster emission 
probability (P n ) for elements with Z > 2 emission from the source 
129 Xe in the l-LGM calculation. The lines show the expected val- 
ues in different temperature with the Poisson assumption according 
to Eq.(2). 



nomenon [Ty,[Tj, we plot ln(n) versus 1/T. Fig. Ogives a 
family of these plots for the disassembly of 129 Xe within 
the framework of l-LGM (left panels) and I-CMD with 
Coulomb interaction (right panels). Again, as Fig.^ the 
clusters are classified according to their masses (upper 
panels), the light isotopes (the middle panels) and the 
charge numbers (the lower panels). For all the panels, 
the obtained Arrhenius plots are linear for the lower T 
side, and their slopes generally increase with increasing 
A or Z value. At these temperatures, one can anticipate 
one large fragment surrounded by many small clusters. 
However, contrary tendency reveals in the high T side 
where ln{n) increases with 1/T, i.e. decreases with in- 
creasing T. In this case, nuclear Arrhenius plots of (n) 
with 1/T are invalid but the Poissonian reducibility still 
remains (see Fig. |2J. This behavior of (n) at higher T 
is related to the branch of the fall of the multiplicity of 
IMF (Njmf) with T where the disassembling system is 
in vaporization [4^, 14a. |44| and hence only the lightest 
clusters are dominated and the liquid residue is vanish- 
ing with increasing T. Above the temperature at which 
the Arrenhius-type plot begins to deviate from the lin- 
earity, the fragments are no longer in coexistence with 
the big liquid drop. In this case, the supply of droplets 



to the vapor has been exhausted and the system is only 
in a single (gas) phase. In other words, the Arrenhius 
law looks valid only when the disassembling system is 
in coexistence phase of liquid and vapor. Recognizing 
this phenomenon, in the following sections we only focus 
on the branch of lower temperature (i.e. in phase coex- 
istence) where the thermal scaling exists to discuss the 
Arrhenius law and their slopes. 



C. "Emission Barriers" 

1. Model dependence 

From Fig. [21 the slope parameter can be directly ex- 
tracted in the lower T side as a function of Z or A. 
In Ref. [I3 Moretto et al. has interpreted these slope 
parameters as "emission barriers" of specific individual 
fragments. This "emission barrier" is usually associated 
with pictures of sequent emission from a compound-like 
nucleus. Fig. Ogives the "emission barrier" of individual 
fragments with different A, ISO and Z in the frame- 
work of l-LGM, I-CMD with/without Coulomb interac- 
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FIG. 3: Arrcnhius-type plot: the average yield per event of differ- 
ent clusters classified with A (top), ISO (middle) and Z (bottom) 
as a function of 1/T. The left panels show the results with I-LGM 
calculation and the right present the results with I-CMD including 
Coulomb force. The solid lines are fits to the calculations using a 
Boltzman factor for (rii). The symbols are illustrated on the figure. 



tion. The error bar in the figure represents the error in 
the extraction of the slope parameter. The first indica- 
tion from this figure is that the " emission barrier" in the 
I-LGM case is the nearly same as the I-CMD case without 
Coulomb force, which supports that I-LGM is equivalent 
to I-CMD without Coulomb interaction rather well when 
the nuclear potential parameter is suitably chosen, but I- 
LGM is a quick model to analyze the behavior of nuclear 
dissociations. The inclusion of long-range Coulomb in- 
teraction makes the "emission barrier" of individual frag- 
ments much lower since the repulsion of Coulomb force 
reduces the attractive role of potential and hence make 
clusters escape easily. The second indication is that the 
"emission barrier" increases with A (Z) at low A (Z) val- 
ues and tend to be saturated at high A (Z) ones. Similar 
experimental trends have been observed for individual 
fragments with different Z in Ref . |15l l45l| or different A 
in Ref. (TT] ■ To make a comparison with the experimental 
data, we plot the data for low energy conditional barriers 
for 82 Kr + 12 C ^ 94 Mo in the same Fig. [J. Since 
the unknown quantitative relation of the temperature be- 
tween the present model and the data, we normalize the 



FIG. 4: The "emission barrier" extracted from the Arrhenius-typc 
plots as a function of cluster mass (top), isotopic mass (middle) 
or cluster charge (bottom) in the cases of I-LGM (solid squares), 
I-CMD without Coulomb (open circles) and with Coulomb (solid 
circles). The unit of the "emission barrier" is MeV throughout this 
paper. The solid lines are fits with the Eq. (8) or (9), and the 
dot-dashed lines represent the fits with the Eq. (10) or (11). The 
open up-triangles and down-triangles are low energy conditional 
barriers for s2 Kr + 12 C — > 94 Mo from Ref. 45j normalized to Z = 
6 for I-LGM case and I-CMD case, respectively. 



above conditional barriers to Z = 6 for I-LGM case and 
I-CMD case, respectively, as done in Ref. [l5j . Overall 
good agreements between the calculation and data are 
obtained. We noticed that the middle panel of Fig. 01 
shows that insensitivity of " emission barrier" of ISO on 
A in the fixed atomic number Z , which indicates that the 
Z dependence of barrier is perhaps more intrinsic and the 
A dependence is basically due to the average effect over 
the species with the same A but different Z. 



2. Source size and Coulomb interaction dependence 

On the origin of these barrier, the surface energy and 
Coulomb energy would play the roles. If the cluster emis- 
sion is mainly controlled by its surface energy, it would 
suggest barriers proportional to Z 2 / 3 (A 2//3 ). In the case 
of I-LGM and I-CMD without Coulomb, we can try to fit 



the barrier for the particles with different mass number 
by 



Bcoui.off = ci x A 2/3 , (8) 

or for the particles with different charge number by 

Bcoui.off = ci x ((A/Z) fit * Zi) 2 ' 3 , (9) 

where (A/Z) fa is a fit coefficient of A/Z for emitted par- 



ticles, and Ai (Zi) is the mass (charge) of particle. c\ is 
the fit constant for surface energy term. The solid line 
in the Fig.Et is a function of Eq.(8) with c\ = 8.469 and 
the solid line in the Fig.^Ji is a function of Eq.(9) with c\ 
= 8.469 and (A/Z) fa = 1.866. These excellent fits imply 
that the surface energy play a major role in controlling 
the cluster emission when the long range Coulomb force 
is not considered. However for the cluster emission with 
the Coulomb field, we can assumed that the barrier is 
mainly constituted by the surface energy term and an 
additional Coulomb term as 



R _ ,2/3 _ 1.44 x Aj/ (A/Z) fit x Z res 

COUL ° n ~ C2X * r Coul (Ai^ + ((A/Z) fit * Z res )V3) 



(10) 
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for the particles classified with different mass number, or 
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Bcoul.oa = c 2 X ((A/Z) f it * Zi) 2 / 3 - 



1.44 x Zi x Z r 



r C oui((Zi * (A/Z) flt )^ + (Z res * (A/Z)j lt y/z) 
I 



(11) 



for the particles classified with different charge number, 
where c 2 is a fit constant for surface term and rcoui is 
chosen to be 1.22 fm. Z res is a fitted average charge 
number of the residue. (A/Z) fa is chosen to be 1.866, as 
taken from the fits for I-LGM. The overall fits for A and 
Z dependent barrier in the case of I-CMD with Coulomb 
force give c 2 = 12.921 and Z res ~ 41 with the dot-dashed 
line in Fig. 0J, and 4c. The excellent fit supports that 
the Coulomb energy plays another important role in the 
cluster emission. 

In the case of I-LGM and I-CMD without Coulomb, 
one would expect the barrier for each Z (A) to be nearly 
independent of the system studied if only the surface en- 
ergy is substantial to the "emission barrier". The left 
panels of the Fig. show the results for Ba, Bjso and 
Bz for three different systems in the I-LGM case. The 
same "freeze-out" density of 0.38po and the same N/Z 



is chosen for the systems of A s 



80 and A s 



274. 



Actually, it appears to have no obvious dependence of 
" emission barrier" on source size as expected for the role 
of surface energy. The solid line in the figure is the same 
as in Fig. 0J However, when the long-range Coulomb in- 
teraction is considered, the "emission barrier" reveals a 
source size dependence. The right panels of Fig. \5\ give 
the "emission barriers" Ba, Bjso and Bz in the case 
of I-CMD with Coulomb force. It looks that the barrier 
increase with the decreasing of charge of system, which 
can be explained with the Eq. (10) and (11) where the 



decreasing of the residue Z res will result in the decreas- 
ing of the Coulomb barrier and hence the increasing of 
the "emission barrier". The lines represent the fits with 
the Eq. (10) and (11) for three different mass systems. 



3. "Freeze-out" density dependence 

In the above calculations, the "freeze-out" density of 
systems is fixed at ~ 0.38po- Considering the "freeze- 
out" density is an important debating variable in the lat- 
ter stage of heavy ion collisions, here we will discuss the 
possible influence of "freeze-out" density on the "emis- 
sion barrier" of clusters. The calculations at the " freeze- 
out" density of 0.177po and 0.597po f° r 129 Xe, corre- 
sponding to 9 3 and 6 3 cubic lattices respectively, are sup- 
plemented to compare. Fig. gives the results of Ba, 
Biso an d Bz at different density. It looks that there are 
no obvious "freeze-out" density dependence in the both 
cases of I-LGM and I-CMD. This is also consistent with 
that assumption that the surface energy is the dominant 
role in controlling the cluster emission. 



IV. SUMMARY AND OUTLOOK 

In conclusion, the Poissonian reducibility and thermal 
scaling of the emitted clusters is explored in the lattice 
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FIG. 5: The source size dependence of the emission barriers for the 
different clusters classified with mass (top), isotopic mass (middle) 
or cluster charge (bottom) from in the cases of l-LGM (left panels), 
I-CMD with Coulomb (right panels). The lines in the left panel are 
fits with the Eq. (8) or (9), and the solid, dot-dashed and dotted 
line in the right panel represents the fits to the "emission barrier" 
of A sys = 80, 129 and 274, respectively, with the Eq. (10) or (11). 



FIG. 6: The "emission barrier" of 129 Xe for the different clusters 
classified with mass (top), isotopic mass (middle) or cluster charge 
(bottom) at the different "freeze-out" density in the cases of l-LGM 
(left panels), I-CMD with Coulomb (right panels). The lines are 
fits with the Eq. (8) or (9) in the left panels and with the Eq. (10) 
or (11 ) in the right panels. 



gas model and molecular dynamical model of Das Gupta 
and Pan. It indicates of a statistical nature of such 
cluster emission. The Poissonian reducibility illustrates 
that the clusters are produced independently each other 
and stochastic. But near the transition temperature at 
which the system leaves the liquid vapor phase coexis- 
tence, the large fluctuation breaks down the Poissonian 
reducibility. The thermal scaling exists while the disas- 
sembling system is at coexistence phase. In this case, one 
can anticipate one large fragment (liquid) surrounded by 
many small clusters (vapor). The calculations are qual- 
itatively consistent with the recent experimental obser- 
vation of Poissonian reducibility and thermal scaling by 
Morctto/Wozniak's group even though the system stud- 
ied is different and the temperature was supposed to be 
proportional to the total transverse energy in their exper- 
iments. This also supports that the lattice gas model and 
classical molecular dynamics is a useful tool to simulate 
the nuclear disassembly. 

Further, based on the Arrhenius law in the liquid 
phase, we extracted the "emission barrier" for clusters 



with the different mass, light isotope mass and charge 
number. The good agreement with the experimental data 
of the low energy conditional barrier is obtained if the 
data was normalization to Z = 6 for the calculation is 
used. Also, the systematic fits with the formula embod- 
ied with the surface energy and Coulomb interaction were 
performed and the overall good fits were reached. The 
results suggest that the cluster emission is mainly con- 
trolled by both the surface energy and the Coulomb in- 
teraction. In the framework of the lattice gas model and 
molecular dynamics model without the Coulomb interac- 
tion, the "emission barrier" relies on the cluster charge 
with the Z 2 / 3 (A 2 / 3 ) law and it does not depend on the 
source size and " freeze-out" density, which indicates that 
the surface energy play a dominant role to control the 
cluster emission. However, in the framework of molecu- 
lar dynamics model with the Coulomb force, the "emis- 
sion barrier" will decrease strongly according to the Eq. 
(10) and (11) and it decreases with the increasing of the 
source size, illustrating that the Coulomb interaction also 
play another weighty role to control the cluster emission. 

Due to the present models are basically statistical, even 
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for the present molecular dynamics model since it uses 
I-LGM predictions as the initial conditions for I-CMD. 
In this case the I-CMD results can not be considered as 
independent ones. Hence, the dynamical transport mod- 
els, such as quantum molecular dynamics [46J, Fermion 
molecular dynamics |47j and anti-symmetrizied molecu- 
lar dynamics [48j etc., will be very valuable to explore 
the statistical and/or dynamical features of cluster emis- 
sion since these dynamical models include some essential 
ingredients of nuclei and their collisions, such as Fermi 
motion, Pauli blocking, mean field effect and collision 
term etc. They can produce clusters by stochastic motion 
without the initial thermal equilibrium assumption and 
"clustering". In fact, some studies have been carried on 
for investigating the binomial scaling of cluster emission 
in terms of an elementary binary decay mechanism us- 
ing a simultaneous statistical multifragmentation model 
and a molecular dynamics model, such as Ref. [4!|. A 
very different feature has been found in that study. In 
the same sprit, it is also meaningful to use a molecu- 
lar dynamics model followed by a statistical decay model 
|50j to explore the Poissionian reducibility and thermal 
scaling. The work along this line is also in progress. 

Finally we would like to emphasize that the present 
work is based on the canonical model calculation, i.e., 
the temperature is fixed in the Monte Cairo simulations, 
which is appropriate for a system in contact with a large 
heat bath. However, a microcanonical treatment to the 
multiplicity fluctuation should be also very interesting 
since the energy conservation will probably play a dif- 
ferent role in fragment fluctuation. Recently, Pratt and 
Das Gupta employed a fragment gas model to discuss the 
importance of energy conservation (canonical vs micro- 
canonical) in affecting these fluctuation [51|. By investi- 
gating a correlation coefficient which was defined as the 
variance of multiplicity distribution relative to its mean 



in such a way that it is positive or negative for super- 
or sub-Poissonian distributions, they found that the cor- 
relation coefficient developed a singularity at the critical 
point when plotted against the temperature in canoni- 
cal calculations; whereas in microcanonical calculations 
the correlation coefficient reached a gentle maximum for 
excitation energies in the fragmentation region. The en- 
ergy conservation (i.e. microcanonical case) results in 
decreasing of multiplicity fluctuation and the multiplicity 
distribution always returns to be sub-Possionian. How- 
ever, even though the singular behavior of the correlation 
coefficient was muted in a microcanonical treatment, the 
behavior of the correlation coefficient as a function of the 
excitation energy was unique to the process of statistical 
fragmentation. Of course, this fluctuation phenomenon 
in different statistical ensembles could be model depen- 
dent. For instance, in a recent paper |52(, Das et. at 
showed that microcanonical calculation of LGM has no 
serious departures from canonical results, eg., the caloric 
curve and fluctuation of IMF distribution are very similar 
between canonical and microcanonical treatment. De- 
tailed analysis and comparisons of the reducibility and 
thermal scaling with the microcanonical LGM are cer- 
tainly welcome. 
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